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Abstract 

This paper compares three different types of "onset of chaos" in the logistic and generahzed 
logistic map: the Feigenbaum attractor at the end of the period doubling bifurcations; the tangent 
bifurcation at the border of the period three window; the transition to chaos in the generalized 
logistic with inflection 1/2 = fJ-Xn ), in which the main bifurcation cascade, as well as the 

bifurcations generated by the periodic windows in the chaotic region, collapse in a single point. The 
occupation number and the Tsallis entropy are studied. The different regimes of convergence to 
the attractor, starting from two kinds of far-from-equilibrium initial conditions, are distinguished 
by the presence or absence of log-log oscillations, by different power-law scalings and by a gap in 
the saturation levels. We show that the escort distribution implicit in the Tsallis entropy may tune 
the log-log oscillations or the crossover times. 
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I. INTRODUCTION 



The transition from regular to chaotic regime presents characteristics similar to phase 
transitions in statistical thermodynamics and an adequate statistical thermo-dynamical for- 
malism has been developed for chaotic systems (see ref. for an introduction). The 
transition point, or border of chaos, is characterized by a null Lyapunov exponent. The 
interesting features of this transition were first illustrated by the pioneering works of Feigen- 
baum on the logistic map attractor at the infinite bifurcation point |^ , immediately followed 
by a series of theoretical works P|, Q], [^. In these investigations a power-law behavior of 
the sensitivity function of the logistic map at the threshold of chaos was formulated and 
;he values for the power-law exponent were analytically calculated. More recently Tsallis 



a 



introduced a formalism of non-extensive thermodynamics that allows to pass easily from 
the chaotic case to the null Lyapunov exponent case, recovering the Boltzmann Gibbs (BG) 
formulation as a limit. This formalism enlightens the connections between chaos and border 
of chaos and, in particular, defines a generalized entropy as the quantity of physical interest. 
For example it has been shown Q that a strong analogy with the Pesin Identity exists 
among this generalized entropy and generalized polynomial sensitivity to initial conditions 
appropriate to describe the sensitivity at the infinite bifurcations point of the logistic map 



Also ot 



groups 



her statistical formalisms emerging within special relativity jl^ and quantum 
seem to describe well the analogies between the threshold of chaos and 



the fully chaotic regime, even if recently the non-extensive formalism has been strongly 
criticized [l3 revealing a still open debate on the subject. At this time a lot of experimental 
results (see IJ] for a review) seem to confirm at least part of the theoretical framework 
of the non-extensive formalism. In this paper we performed various numerical experiments 
devoted to investigate the logistic map at transition points presenting different routes to 
chaos. Thus our study regards different border of chaos regimes with different attractors 
underlying the dynamics. We present a joint analysis of (generalized, coarse grained) entropy 
and occupation number. We investigate the dynamics at large times and at large sampling 
ratios and analyze border of chaos transitions different from the Feigenbaum attractor in 
the logistic map. Finally we examine the presence of oscillations in the convergence to the 
attractor. 
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II. NUMERICAL EXPERIMENTS ON THE FEIGENBAUM ATTRACTOR 



In this paper we describe a set of numerical experiments, performed on the logistic map 
at the onset of chaos, with the purpose of examining the dynamics of a statistical ensemble 
of points, starting from far-from-equilibrium initial conditions. In each experiment the 
available phase space is partitioned in a number Wbox of elementary, equal cells. A set of N 
points is randomly selected in the phase space according to two kinds of initial set-ups, and 
iterated according to the map. In these experiments we used the logistic map in the form 

Xi+i = 1 - /ix^ ; -1 < a; < 1 ; < < 2 . (1) 

We investigated the parameter values n = 2 for the fully chaotic 1.401155189.. at 

the Feigenbaum attractor, jitg = 1.75 at the tangent bifurcation and /ii/2 = when the 

1 /2 

logistic is generalized with inflection 1/2 {xi+i = 1 — fix/ ). The first kind of initial start-up 
is from concentrated initial conditions (i. c). All the N points are chosen, with a uniform 
random distribution, inside a single cell of the partition, itself chosen at random. In the 
second kind of experiment all the N initial points are uniformly and randomly distributed 
in all the available phase space (spread i. c). The i-th cell will contain a fraction Nj/N of 
the total number of points (Yl^r ~ that one can naively define a probability 

of occupation for the i-th cell through pi = Ni/N with the constraint Pi ~ 

one hand we observe, for both kinds of experiments (concentrated and spread i. c.) the 
occupation number in time, namely the number of non empty cells of the partition. On the 
other hand, the other quantity of interest in our discussion is the physical entropy, defined 
as coarsegrained entropy through the probabilities pi. We will use the Tsallis definition of 
entropy iQ], useful to characterize the behavior at the onset of chaos where the Lyapunov 
exponent vanishes: 

T _ Sr'^box <i 

qTS _ ^ 2^1=1 Pi fo\ 

Other definitions of entropies Il2l | would also work for the purposes of this paper. 

The analysis is performed through comparisons of the various behaviors between chaotic 
regime and border of chaos. We use different partitions Wbox and different sampling ratios 
r, defined as r='N /Wbox, that gives an indication of the goodness of sampling. 

The fact that the Lyapunov exponent is zero at the threshold of chaos, allows us to use 
power-laws to characterize the evolution P, U, 3|- The comparison with the chaotic 



case makes natural to pass from an "ex pon ential formalism" to an "extended exponential 
formalism" describing the power-law [tI. llT^. As discussed in 3] , in the Chaotic regime, 
the (averaged) time evolution of the map, starting from concentrated initial conditions, may 
happen in two or three stages: a first "thermalization" stage heavily dependent on the details 
of the initial conditions; a second "linear" stage in which the coarse grained Boltzmann-Gibbs 
(BG) entropy Sbg grows linearly with time and a third "saturation" stage in which Sbg 
reaches its final equilibrium value. 

In the fully chaotic regime, when proper average is made on different experiments of the 
same kind (see ref. 3]), the evolution of the logistic map, starting from concentrated initial 
conditions (CIC), exhibits an initial exponential growth of the occupation number without 
a first thermalization stage. A saturation is reached, depending on the grid size, when Nocc 
equals the fractal support value (Nocc = ^box for /i = 2). The coarse grained version of 
the BG entropy shows a linear initial increase followed by a saturation at a S^Csat level 
i^BGsat — ^^(Wbox)) when starting from CIC. A Pesin-like identity can be observed for the 
coarse-grained BG entropy. Using for this case a non-extensive formalism 0, 12 , , a 

linear growth of Sg in the first evolution stage can be found and a Pesin-like identity can be 
recovered, if the appropriate entropic index is selected . Thus some features of the chaotic 
case can be directly transferred to the threshold of chaos. It is important to note that in 
the chaotic case the evolution happens in only two stages. In this paper we are interested 
in clarify how the final saturation stage is reached at the onset of chaos. To this aim we 
studied the evolution at longer times than the previous works. The results are illustrated in 
FigsQandEl where, differently from the chaotic case, we can observe roughly three stages. 
The first stage is the analogue of the chaotic one: both {Nocc) and (Sg) increase and reach 



a maximum. The duration of this first stage becomes longer and the maximum 



12 



evel higher 



19|, in this 



increasing Wbox, while they do not depend on r. As previously stressed 0, 
first "linear" stage, (Sg) grows linearly only for the proper entropic index {q = 0.36 when 
A* = /^oo)- In the second stage both {Nocc) and (Sg) decrease. (Nocc) follows a power-law with 
superimposed log-log oscillations which amplitude increases with r and does not depend on 
Wbox- For low enough q values, (Sg) reproduces the {Nocc) behavior (when g ~ results 
Sg ~ Nocc "1)? but the amplitude of oscillations decreases with increasing q and vanishes 
when q ~ 0.36 (see figO right frame). In the third stage {Nocc) and (5*^) reach their (Wbox 
dependent) saturation values and the evolution ends. Note that the saturation is reached 
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FIG. 1: (a): time evolution of the average occupation number {Nqcc) for the Logistic map at /ioo- Different 
curves correspond to different Wbox with r—1 (from top to bottom Wbox = 10^, 5 • fO'', 2 • fO'*, fO''). 
(b): Nocc vs. time from spread initial conditions (upper curves) and (Nocc) from CIC (lower curves). Fixed 
Wbox (10^) and different r (from top to bottom: r— 1000,100,10,1 upper curves, r— 100, 10, 1 lower curves). 
The averages are made over 5000 randomly chosen initial configurations. 

at a level that increases with Wbox, but (Nocc) results in a lower final value than the fractal 
support, that is equal to Nocc = 244 (799) for Wbox = 10^ (10^). 

The log-log oscillations of Nocc have been already observed in the logistic map at = fi^o 



22], performing numerical experiments starting from uniform i. c. (see fig^-b). Here, 



for the first time, we can observe this same behavior for (Nocc) starting from CIC (see 
figffl-b). This kind of logarithmic oscillations, superimposed to a power-law, characterize 
a large number of systems exhibiting discrete scale invariance j23j. The time evolution 
Noccit) = t^^-P(ln(t)) (where P(ln(t)) is a periodic function) can be expanded in Fourier 
series and, keeping only the first term, can be written in the form: 



Nocc{t) = t-^ A{1 + Bcos{u\n{t) +</>)) 



(3) 



As can be observed in fig^b (upper curves), starting from uniform i. c, the exponent 6 in- 
creases with r (this dependence has been stressed in j22|) and its limit value {S — 0.800138... 
for r 00) has been determined in We analyzed the results of our experiments from 
concentrated i. c, using for {Nocc)(t) the same time dependence of eq.©, founding the 
power-law exponent 6 ^ 0.54, but we could not establish any dependence on r. Note that 




FIG. 2: Logistic map at the critical point /i = 1.401155189. Tsallis Entropy (Sq) is averaged over Nav = 5000 
random choices of an expansion experiment initial cell and is plotted together its statistical error, (a): 
different point ensembles are relative to increasing Wbox values (from top to bottom Wbox — 5 • 10^,2 • 
10'"', 10^, 5 • 10**, 2 • 10^,10"') while the entropic index q = 0.36 (the value that linearize the first evolution 
stage) and the sampling ratio r = 1 do not vary, (b): fixed Wbox = 10'* and r = 100, the point ensembles 
are relative to different entropic indexes (from top to bottom q — 0.001, 0.1, 0.2445, 0.36, 0.7, 1). 

also starting from uniform i. c. the evolution can be roughly divided into three stages: a first 
"crossover time" stage, in which Nocdf) —const. ; the power-law stage, with superimposed 
log-log oscillations; the final saturation stage, in which Nocc{t) equals the fractal support 
value, differently to the CIC experiments. 

To investigate the origin of the observed behavior (on average) starting from CIC, we 
performed others expansion experiments, without averaging, selecting a special initial cell. 
In figiniwe show the results of an expansion experiment from the VFbox-th cell, that has the 
point X = 1 as right extreme, where the distance of x = 1 from the Feigenbaum attractor 
is infinitesimal. Again there are three evolution stages. The first stage presents, for both 
Nocc{t) and Sq, exactly the same pattern already observed for the sensitivity function starting 
from xq = 1, characterized by large fluctuations 0]. The upper limit values correspond to 
times n = 2'^ — 1 and form a straight line for Noccit) and S'g=o.2445 as already found for 
the sensitivity function |9]. This stage is longer when Wbox increases. In the second stage, 
the upper limit of the large fluctuation pattern decreases, following approximately a power- 
law with superimposed log-log oscillations. Finally a saturation stage is reached, in which 
the evolution shows a regular periodic pattern. Some features of the evolution from this 
special cell remind the averaged behavior. This is not trivial because there are single cells 
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FIG. 3: Time evolution of the Logistic map at the critical point /i = 1.401155189, the maximum values 
are reached at time n = 2'"' — 1 (with k — 1,2, 3, ...) and are stressed by circles, minimum in 2*^ by squares. 
Nocc vs. time obtained by an expansion experiment starting from an initial condition concentrated in the 
extremely right cell (centered in = 1 — l/{2Wi,ox)): sampling ratio r = 1, number of partition cells 
Wbox = 10^ in (a) and Wbox = 10^ in (b). 

which exhibit a completely different time evolution. We check;ed, for instance, the repellor 
X = Xr = 0.5602326... (solution of the equation Xr = I — fioo^r)- Centering the initial cell 
in Xr, we observed a first stage in which Nocc grows exponentially and reaches its maximum 
at Nocc ~ Wbox] then it decreases following a power-law without log-log oscillations and its 
evolution ends at a saturation value equal to the fractal support. Large fluctuations of Nocc 
and Sq do not appear at any time. 

The Feigenbaum-scaling cascade of period doubling is not the only possible scenario. In 
the following we show few comparisons with others routes to chaos. The logistic map at the 
beginning of the period three window {fi = fit = 7/4) exhibit both a vanishing Lyapunov 
exponent (A = 0) and a power-law decreasing sensitivity function (weak insensitivity). The 
power-law exponent and the correspondent Sg entropic index (g = 3/2) have been analyt- 



ically determined in 25|. We performed for the tangent point fi = fit the same relaxation 



experiments described in the previous part for fi = fi^o, reproducing and improving the 
results already obtained in and {2^. Our results are showed in figlH As already stressed 



in [2a|, the evolution starting from CIC begins miming a chaotic behavior: Sbg grows lin- 
early and (Nocc) exponentially; then a maximum is reached. Thereafter (Nocc) presents a 
"crossover" stage, with (Nocc) = const. ~ W^ox, which is longer for larger r. This has been 
already analyzed in |j2J] for uniform i. c. but not for CIC. In the third stage (Nocc) decreases 
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FIG. 4: Logistic map at the tangent point /i = fit- Wbox = 10^, averages performed using Nav=800 random 
choices for the initial cell, (a): Nqcc from a relaxation experiment starting from uniform i. c. (continuous 
lines) and {Nqcc) obtained starting by a CIC (points), are plotted together for different r values (from top 
to bottom r = 1, 10, 100). (b): Tsallis entropy with index g = 3/2 is showed for different r (top to bottom 
r = 1, 10, 100); the three series of points are inside their statistical errors (not showed in figure). 

with a power-law and no log-log oscillations. The evolution ends when (Nocc) reaches the 
three points attractor. The curves relatives to CIC and to uniform i. c. join each other after 
the first stage (fig. ID -a), differently from the (figffl-b). Turning our attention 

to fig. El-b Sq, after the mimed "chaotic" stage, decreases and is linear only for q = 3/2. 
There i= no co^ove. t.n,e fa any sampling ratio and for , ^ 3/2 the dope depends on 
|26| but not on r. 

Finally we examined the transition to chaos in the generalized logistic map when the 
inflection is 1/2. Here the main bifurcation cascade, as well as the bifurcations generated 

n 

by the periodic windows in the chaotic region, collapse in a single point [27]. In such a 
point the map undergoes the transition to chaos and the Lyapunov exponent is zero. The 
experiments performed starting from uniformly spread i. c. show again two features already 
encountered in the previous cases. There exists a "crossover time " (figEl) which is more 
extended for larger r. Thereafter it appears a regime of convergence to the attractor in 
which Nocc follows a power-law with negative exponent and no log-log oscillations. In the 
limits of our numerical experiments such exponent depends neither on the grid nor on the 
sampling ratio. 
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FIG. 5: Evolution from spread i. c. at ^1/2 in the generalized logistic with inflection 1/2. The curves are 
obtained, from top to bottom, using r = 100, 10, 1 for W^ox = 64.000 (a) and W^ox = 128.000 (b). 

III. CONCLUSIONS 

Our numerical experiments show that the long time evolution on the Feigenbaum attrac- 
tor from CIC happens in three stages, instead of the two observed in the chaotic case. The 
saturation is reached only after a power law decreasing stage, with superimposed log-log 
oscillations, observed for high enough r values. 

We also showed that crossover-time appears in the time evolution from uniform i. c, in 
all the cases examined. The crossover-time appears also in the time evolution from CIC of 
the logistic map for fi = fit but not for /i = /ioo- 

Observing Sq in the expansion experiments from CIC for high enough values of the 
entropic index q, log-log oscillations (for /i = /ioo) and crossover-time (for /i = fit) do not 
appear. The escort distribution, implicit in the Tsallis entropy formulation eq.©, selects, 
increasing q, the more populated regions of the phase space. Thus the two observed behaviors 
probably originate from the contribution of the less populated cells. 
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